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ABSTRACT 

A  review  of  recent  results  associated  with  e*actly  solvable 
multidimensional  nonlinear  systems  and  related  duestions  of  direct 
and  inverse  scattering  is  given. 
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In  tnis  lecture  a  review  of  some  recc-nt  results  associated 
exactly  solvable  multidimensional  nonlinear  systems  will  ee  given. 

Tie  motivation  for  much  of  this  work  has  come  via  what  '$  cownonly 
referred  to  as  tne  Inverse  Scattering  Transform  (I.S.T.,  as  a 
reference  See.  fo<  examp le .  1ST  is  3  method  to  solve  vertain 
nonlinear  eauaticns  by  associating  the*  with  appropriate  compatible 
J 1  near  equations,  :ne  of  -men  is  identified  as  a  scattering  proD'.e’” 
ana  tne  otners(s)  serves  to  fix  theVtime  evolution"®^  the  scattering 
data.  — — 

In  one  spatial  dimension  tne  prototype  problem  is  the  (KdV) 
equation 


The  KdV  equation  is  compatible  with 

vxx  ♦  u(x,t)v  *  »v 

vt  *  ^Y'*UI^V  *  (-•*2u)»x 


v  .  *  V. 

xxt  txx 


implies  (I).  Equation  .2)  is  the  Schrodinger 


scattering  proolem,  *  the  eigenvalue  (  ■»*  const,  in  (3)).  The 
solution  of  (1)  on  the  line:  <x<»for  initial  values  u(x,t*0) 
vanishing  sufficiently  rapidly  at  infinity  is  ODtained  by  Studying 
associated  direct  and  inverse  scattering  problem  of  (2)  anq  us’rg  ; j 
to  fix  the  time  evolution  of  the  scattering  data.  It  turns  out  that 
the  inverse  problem  amounts  to  solving  a  matrix  Riemann-Hi Ibert 
boundary  value  problem  (RH8VP)  whose  jump  discontinuity  depends 
explicitly  on  the  scattering  data.  Calling  * «-k2 ,v( x  ,k ) *u ( x , x )e* 1 
the  RHflVP  takes  tne  following  form, 

(u^-w  J(x,t,k)  *  u_(x,t.c(k))  V  (  x  ,  t ,  k  )  on 


u:“l,  |  k  |  — 


where 


V(x.t.k)  -  r(k.t)  e2,lt\  a(k)  «  -k,  :•(  k-.ktli  ,  and  u.  are  the 

limiting  boundary  values  as  ImkiO:  of  meromorphic  functions  in  the 
upper  (♦)  lower  (-)  half  plane.  (4)  may  be  converted  into  a  linear 
integral  equation  by  taking  a  minus  projection  and  the  potential  is 
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reconstructed  via 


j;».t)  *  -  -  rj  /  L (k.x.t.-k)  V(»,:,»jak  (5) 

where  the  contour  is  taxen  aoove  all  Doles  of  r(k,t).  of  -men  thiji-u 

IS  at  "lost  *  ,,r'tC  ru‘'ipfrr ,  k  «  )k  ,  v  '0  J  s  '• M.  Thy 

J  J  J 

scattering  oa;a.  fr  reflection  coefficient,  n.k.t)  evolves  simoly 
in  time  ^ 

rU.t)  *  n(k.O)  e8lk  :  16) 

The  above  sene**  may  oe  extended  so  as  to  solve  a  surprisingly 
large  number  of  interesting  nonlinear  evolution  equations.  There  are 
two  scattering  proplems  pf  particular  interest  in  one  dimension: 

(i)  Scalar  scattering  problems: 

*  r  u  { * )  dn~^v 

57  j-2  J  ST’*7  *  kv‘  ‘ 

v( x ,k) ,  Uj  t  ( 

(n)  first  oraer  systems  -  generalized  akns 
—  *1  k  J  v  *  Q  v 

OX 

v ( x , k ) . q( x )  c  tNxN,  J  *  diag  (J^.*--Jn) 
j’*JJ  .ifj 

qn  .  0. 


Via  an  appropriate  transformation  the  inverse  problem  associate  with 
(i).  (ii)  can  be  expressed  as  a  matrix  RHBVP  of  the  form  (4).  The 
potentials  u,,q  can  be  shown  to  satisfy  nonlinear  evlution  equations 

J 

by  appending  to  (t),  (u)  suitable  linear  time  evolution  eouations. 

One  then  finds  that  the  scattering  data  V(x.t.k)  evolves  simply  in 
time.  Well  known  solvable  nonlinear  equations  include  the  Boussmesq. 
modified  KdV,  sine-Gordon,  nonlinear  Schrodinger,  and  three  wave  in¬ 
teraction  eouations.  The  reader  may  wish  to  consult  for  example^4"*^ 
for  a  detailed  discussion  of  some  of  this  material. 

It  is  most  significant  that  these  concepts  can  pe  generalized  to 
2  spatial  plus  one  time  dimension.  Here  the  prototype  euqation  is 
the  lUdomtsev-Petv iashv il i  (k-P)  equation: 
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which  is  the  compatibility  equation  Between  the  Allowing  linear  pr;;. 


v«x  *  u(*.y.t)v  ‘  0 


V.  *  4v,,,  *  6uv>  *  3( i  %dx  '  )  v  »  >V  *  0  (?) 

(Y  •  const.)-  we  Shall  consider  the  question  of  solving  (7)  for 

2  2  2 

u(x,y,0)  decaying  sufficiently  rapidly  in  the  plane  r  »  ♦  y*  .  «. 

Physically  speaking,  Both  cases  •  -1  (KPI)  *  ♦!  (KPII)  are  of 

interest.  Whereas  KPI  can  Be  related  to  a  RHBVP  of  a  certain  type 

(nonlocal;  see  ref.'*')  KP11  turns  out  to  require  new  ideas.  Letting 

2 

,  i  kx  ♦  k  y/O 

v  *  . ( * ,y ,k)e  1 

a  m  o  Q  *  iOp  .-g  t  0.  Then  there  exist  functions  u  bounded  for  all 
x,y  satisfying  .  -  1  as  |k  |  -  However  Such  a  function  turns  out 

to  be  nowhere  analytic  in  k,  rather  it  depends  nontrivially  on  Both 
the  real  and  imaginery  parts  of  k  «  (kR  *  ikj).  -  *  u(x,y,kp,k.) . 

In  fact  v  satisfies  a  general ization  of  a  PHBVP  •  namely  a 
I  ( D8AR)  problem  where  _  satisfies, 

«  uU.y.Cg.kj)  V ( x ,y , kR ,  k , ) 

where  —  ■  i(Tr-  ♦  i  =-§— )  and  V  has  the  structure 
**  "  1  ifi(x.y.kR,kj ,£q) 


V(x,y,kp,ky)  * 


sgn(kg)e 


T( kR,k j ) 


X  X. 

i ( x ,y , kp , x . , ir, )  1  (x  *  2 y  *  kp)  *  *2(x  *  2y  r^lkg 

Jc  *  *  7^kr  ko  ’  kR  *  ^ki  111 

(11)  may  be  converted  into  a  linear  integral  equation  by  employing 
the  generalized  Cauchy  formula.  T ( k R , k j )  is  viewed  as  the  ("nonphysi¬ 
cal"  data,  i.e.  inverse  scattering  data:  i.e.  inverse  data)  and  the 
potential  is  reconstructed  via 


F//  j  JJ  * 


• "  V*«_- 


U  {  1  ,y  )  *  *  r—  j  j  U\X,y,£Q,kj)V^X,y,k^,k|)dk£dkj.  (12, 

For  k-P  tr>e  evolut'on  of  the  data  obeys  (y  *  i,k^  ,n  (9)) 

~  *  (8'Y. ;(►*»-  -  -k t  •  3k^)T  'll 

■  t  -  u  0 

-here  k0  ’  kR  *  *“  •  k  *  kR  *  ,kI- 

Similar  ideas  apply  to  higher  order  scal*r  problems 
„n  n  .n-j 

(Hi)  :°JL  .  LJL  .  :  u.(jl)2 — IX  .  0 

**  3xn  j*2  J  3.n*J 

where:  v.  Uj  f  t  and  to  First  order  systems 
(iv)  ;  *  ~:rf  *  0(*.y)v  *  0 

cy  9a 

where:  v,qcCNxN,  Jnliag!  J1 , . . .  ,JN) ,  j’  ^  JJ  .  i  P  j  with  q' ’  *  0. 
Interested  readers  many  consult  reference4* For  associated  details. 

The  notion  of  a  extends  to  higner  dimensional  scattering  and  in¬ 
verse  scattering  problems.  However  as  «e  shall  mention,  oespite  the 
Fact  that  the  .inverse  scattering  problem  is  essentially  tractable 
there  does  not  appear  to  be  any  local  nonlnear  evolution  equations  m 
dimensions  greater  than  2*1  associated  wit*  *hj1  tidimensional  gener¬ 
alizations  of  (iii)  or  (iv). 

Our  prototype  scattering  problem  will  be 

ov  *Av*u(x,y)v*0 
n  .2 

1  *  :  2_,  X  t  Rn.  y  e  R  .  (14) 

1  «1  axf 


letting 


v  *  _(x,y,k)e 


ik-x  *  k  y/c 


k  •  kR  *  ikp  k  £  1 

k-x  .  :  kjXj.  O  •  cfi  *  iCj. 

Then  there  exist  functions  u  bounded  for  all  x,  y  satisfying  .-1.  as 

|k.  |  -  j  •  1 . n.  when  c -  t  0  u  turns  Out  to  be  nonanalytic 

J  " 

each  of  the  variables  k,  i.e.  u  *  u(*,y,kfi  ,...kR  .  k.  . k;  )  and  . 

1  n  1  n 


/  .*  • 
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satisfies  a  ;•  problem  linear  in  ...  m  each  of  the  variate  k  ,  -.a. 

J 

satisfies  an  equation  of  the  form, 

*  T^(.  ) ;  j  *  1 . n  (  1; ; 

J 

where  T.(u)  is  an  appropriate  linear  integral  operator  which  depends 
J  •  " 

only  on  one  scalar  scattering  function  T  ;  T\  ■  T^[T],  T*T(kp,kj,C) 

£  being  (n-1)  integration  parameters  in  the  nonlocal  operator  T . .  The 
inverse  problem  is  redundant,  i.e.  we  are  given  T(kR,kj,£)(3n-i)  para¬ 
meters)  and  we  must  reconstruct  a  local  potential  u(x,y)(n+l  para¬ 
meters).  A  serious  issue  is  how  to  characterize  admissible  inverse 
data  T,  i.e.  data  that  really  arises  from  a  local  potential  (small 
generic  changes  in  T(kR,kj,{)  cannot  be  expected  to  arise  from  a 
local  potential  u(x,y)).  Insight  into  this  question  is  obtained  by 
requiring  (1  t  j).  The  form  of  this 

constraint  is  given  by 


*,j(T)  •  Ni J[T]  (16) 

where  i  is  a  linear  operator  and  N. .  a  nonlinear  (quadratic)  non- 

J  'J  c,|,] 

local  operator.  Details  can  be  found  in  ’  .  Equation  (16)  can  oe 

integrated  and  this  integrated  version  may  be  used  to  reconstruct 

u(x,y)  as  well  as  give  a  characterization  for  admissible  scattering 

data:  T(kRtkj,0.  However  (16)  also  indicates  why  simple  local 

nonlinear  evolution  equation  have  not  been  associated  with  equation 

(8).  Namely  in  the  previous  lower  dimensional  (2*1  and  1*1)  problems 

the  time  evolution  of  the  scattering  data  obeyed  a  particularly  simple 

equation,  (e.g.  4^  ■ui(kR,kj)T.  However  in  this  case  such  a  simple  flow 

will  not  be  maintained  -  due  to  the  nonlinear  constraint  (16). 


These  ideas  can  be  generalized  to  first  order  systems: 

i  n  , 

(v)  : —  ♦  cl  J-  i—  *  qv 

J-I  J  “j 


kM. 
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.ijn  similar  results  oota 1  ned6a ’ 6  y  Again  tne  scattering  cata  satis¬ 
fies  a  nonlinear  constraint.  !n  general,  tnere  is  no  compatible  local 
nonlinear  evolution  equation  associated  with  (v).  However  when  cer¬ 
tain  restrictions  are  out  on  tnen  tne  constraint  equation  oecomes 
linear  ano  tne  so-called  n  wave  interaction  equations  are  compatible 
witn  the  system  ( v ) .  Nachman  and  Adowi  tz^a  4  showed  tnat  at  most.  ; re 
system  would  De  3*1  dimensional,  and  Fokas^0^  showed  tnat  indeed  tne 
system  is  reducible  to  2*1  dimensions  by  a  transformation  of  indeoend- 
ent  variables  (cnaracteristic  variables). 

Beals  and  Coifman  have  given  an  alternative  but  similar  formgla- 
7a  hi 

tion  a,0J  m  the  scalar  case. 

There  is  an  n-oimens lonai  problem  which  also  fits  within 
the  framework  Of  1ST:  The  so-called  generalized  wave  and  generalized 
sme-Gordon  euqation  (GWE  and  GSGE).  These  equations  arise  in  the 
context  of  differential  geometry  and  serve  to  extend  the  classical  re- 
suits  of  Backlund  for  the  sine-Gordon  equation  to  n-dimensions 
The  n-dimensional  Backlund  tranformati on  is  given  by: 

dX  *  XA:x  «  A  -  XB.  ( -:1 


..  -  5*  - 

dx  *  _  - —  dx , , 

J-l  J 

Aij  *  6i<z)VV 


l  “•l: 


l  -ii 


5  2  -  t 1  dx  .  -  - —  -r - dx  ,  1  <  i »  J  <  n ,  (18) 

ij  an  ,xi  j  Au  1  " 

and  a  *  (a.  :  c  Fnxn.  Equations  (17-18)  reduce  to  tne  Backlund  trans¬ 
formation  for  the  generalized  sine-Gordon  equation  (GSGE)  when 

£,(z)  *  (z2  *  (24n  -  1 ) )/2z ,  (19) 

and  for  the  generalized  wave  equation  (GWE)  when 

S,(z)  *  -(l-z2)/2z  i  x(z).  (20) 

The  compatibility  condition  required  for  the  existence  of  solu¬ 
tions  to  these  Backlund  transformations  results  in  a  system  of  second-. 
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oroer  partial  differential  equations  for  an  ortnpconal  n  «  r.  -natr-, 

a  1  •  a 7 j i  m  (17)  »mcn  is  a  function  of  n  independent  variates 

a  =  a(x,,x, . *  ).  The  equation  nas  the  form 

1  l  n 


i-f J-i 

i  V  li  i 


__  /  J__  °1 1 

3\)y 

1  Sali  ‘’“it 


k  f~'’i  «ik 


=  ..  a , ,  a 


•-an‘ir 


i  5iH 


_ 1_  ^lk 

alk4lj  3*k  3*j 


i ,  J  ,  k  distinct. 


!iii  ,°hi  ill!  ,  ,  i. 

;*k  *u  s>i 

where  c  ■  1  for  the  GSGE  and  c  *  0  for  the  Gw£. 

We  observe  that  when  n  «  2  and  c  1  1  (GSGE),  the  orthogonal 
matrix  a  *  { a  ■  s }  given  by 


cos  2  u 


Sin  2  u 


CCS  *  u 
c 


for  the  function  u  *  u(*,t)  reduces  the  GSGE  to  the  classical  sine- 
Gordon  equation  {<  *  -1), 

utt  '  u*x*<sin  u  ‘  °-  (23) 

On  the  other  hand  wnen  n  *  2  and  <  1  0,  then  with  (22)  the  GWE  reduces 
to  the  wave  eauation  (23).  When  n  »  3  the  general ization  of  the 
wave  equations  discussed  here  is  nonlinear. 

The  Mckl und  transformations  (17)  described  above  are  in  fact 

matrix  Rlccatl  equations.  Linearizations  of  such  a  system  can  be 

9] 

performed  in  a  stnaghtforward  manner  (see  for  example  J)-  Intro¬ 
ducing  the  transformation 

X  *  UV*1,  (24) 

Where  U,  V  and  n  X  n  matrix  functions  of  The  following  linear 


^  ,.Vv 


system  is  deduced: 


with  the  components  cf  A,  B  given  by  (18).  Compatibility  ensures  t •'a * 
the  orthogonal  matrix  a  *  ;a-ji  satisfies  cue  GSGE  with  (15)  ana  G»E 


* 1 1 m  (20).  Al ternat i vely ,  if  we  call 


(:)•*• 


the  following  linear  system  of  2n  o.d.e.'s  are  obtained: 

.■  ^ *  X  A  •  Ui  ♦  C-lfc, 

«*j  J  J 

where  A  ,  C  are  2n  x  2n  matrices  with  the  block  structure 
J  j 


Here  2j,  are  n  x  n  matrices  having  the  following  structure: 


i.  *  (-  -  1  )e,a  .  *  a  , 

J  *  »  J  J 


4J  *  4eJ 

where  e-  *  {e .  > - .  is  the  unit  matrix 
J  J  f 


{ej}  ik 


i  *  k  *  j , 


otherwi se , 


and  in  component  form  takes  the  form 

'■At  • 


4  I,  .  _J 

lj  ‘  iJ;»u  «| 


In  (28)  a  is  the  orthogonal  matrix  Rn  -  S0(n)  associated  with  the  GhE 

when  6  ■  X  and  with  the  GSGE  when  6  •  |(z  ♦  1/z),  *«  ^(z  -  1/z).  and 

y,  is  the  matrix  (30):  F  -  *  (F),  y,  ♦  >  «  0.  Eouations  (21)  arise 
J  n  n  j  j 

as  the  compatibility  condition  associated  with  (26).  More  explicit*, 
for  the  GWE  the  scattering  problem  takes  the  form  [*  «  *(*.*)] 

j  £  »  x  a  a*  ♦  C  ii*  (31) 

«*j  J  J 


.  T , <T.  r.  «.-« »r.  V,  K-.  T. 


V 


a 


and  C;  given  py  (27,30) . 

W 

for  tne  GSGE  tne  scattering  proplem  for 


*  c(; 


i «.:)  taxes  t-e 


n-W, 


**  Cjfe 


Hz),  *(z] 


given  aoove,  or  equivalently 


iTT’i  Vf8j ' 


♦  CjO. 


u  *  diag(*l,  -1, 


8l 

In  J  it  is  shown  how  these  linear  problems  may  be  viewed  as  a  direct 
arc  inverse  scattering  problem  for  the  GWE  and  GSGE.  Namely  tne 
direct  and  inverse  problem  may  be  solved  for  matrix  potentials,  de¬ 
pending  on  the  orthogonal  matrix  a,  tending  to  the  identity  sufficient¬ 
ly  fast  in  certain  "generic"  directions.  It  should  be  noted  that 
solving  the  n-dimens iona 1  GWE  and  GSGE  reduces  to  the  study  of  the 
scattering  and  inverse  scattering  associated  with  a  coupled  system  c* 
n  one-dimensional  o.d.e.'S.  This  is  in  marked  contrast  to  other 
attempts  described  earlier  to  isolate  solvable  (local)  multidimensional 
nonlinear  evolution  equation  which  are  compatibility  conditions  of  two 
Lax-type  operators,  e.g., 


where  lisa  partial  differentia!  operator  with  tne  variable  t  enter¬ 
ing  only  parametrical ly.  Although  as  we  have  seen  nonlinear  evolution 


1 
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equations  in  tnree  independent  variables  can  De  associatec  wtn  s-c" 
tax  pams  (e.g.  tne  k-P ,  Davey-Stewartson ,  tnree  wave  interaction 
equations,  etc.)  little  progress  via  tins  route  nas  seen  made  in  mc'-e 
tnan  tn>-ee  dimensions.  As  discussed  earlier  one  nas  to  overcome  a 
serious  constraint  mnerent  in  tne  scatter i ng/ 1 nverse  scattering 
theory  for  mgner  dimensional  partial  differential  operators  m 
order  to  be  able  to  isolate  associated  solvable  nonlinear  equations, 
i.e.  the  scattering  data  generally  satisfies  a  nonlinear  equation 
(e.g.  (16)).  Tne  analysis  associated  with  the  GWE  and  GSGE  avoids 
these  difficulties  since  the  GWE  and  GSGE  problems  are  simply  a 
compatible  set  of  nonlinear  one-dimensional  o.d.e.'s.  The  results 
in  [8]  demonstrate  that  the  initial  value  problem  is  posed  with 
given  data  along  lines  and  not  on  (n-1)  dimensional  manifolds. 

Similar  ideas  apply  to  certain  n-dimensional  extensions  of  the  so- 
called  anti -sel f-dual  Yang-Mills  equations  (SDYM).  In^  it  is 
shown  that  these  multi-dimensional  nonlinear  equations  are  associated 
with  compatible  two-dimensional  linear  systems.  Broad  classes  of 
solutions  may  be  calculated  by  tne  c  method.  Since  the  overall  com¬ 
patible  linear  systems  are  coupled  two-dimensional  equations,  the 
scattering  data  does  not  satisfy  the  nonlinear  constraint  discussed 
ear  1 i er. 

Finally  we  remark  that  there  is  a  class  of  nonlocal  equations 
which  can  be  reduced  to  exactly  solvable  equations.  In  the  context  of 
multidimensional  nonlinear  equations  perhaps  the  most  interesting 
example  is 

(ut  *  Uxxx  *  2(uH:u)x}x  =  -3c%y’  ]£! 


u  W  v-~. 


[H,u)ix,y,2.t)  *  J  *  uLLJi^)  dE 


and  i  denotes  the  Cauchy  principal  value  integral.  (38)  is  reduced 


to  the  K-P  euqation 


( w  *  w  -  1  (w2 )  )  *  -  3c2w 

t  XXX  x'x  yy 


via  the  transformation 


U'lH^U. 


Details  and  other  examples  are  given  in  [10]. 
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